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ABSTRACT 

We present a statistical analysis of a sample of 20 strong lensing clusters drawn from 
the Local Cluster Substructure Survey (LoCuSS), based on high resolution Hubble 
Space Telescope (HST) imaging of the cluster cores and follow-up spectroscopic ob- 
servations using the Keck-I telescope. We use detailed parameterized models of the 
mass distribution in the cluster cores, to measure the total cluster mass and fraction 
of that mass associated with substructures within R ^ 250kpc. These measurements 
rS \ are compared with the distribution of baryons in the cores, as traced by the old 

■ stellar populations and the X-ray emitting intracluster medium. Our main results in- 

clude: (i) the distribution of Einstein radii is log-normal, with a peak and la width 
of (\og 1Q 0E(z = 2)) = 1.16 ± 0.28; (ii) we detect an X-ray/lensing mass discrepancy 
of (Msl/M x ) = 1.3 at 3cr significance - clusters with larger substructure fractions 
displaying greater mass discrepancies, and thus greater departures from hydrostatic 
equilibrium; (iii) cluster substructure fraction is also correlated with the slope of the 
gas density profile on small scales, implying a connection between cluster-cluster merg- 
ers and gas cooling. Overall our results are consistent with the view that cluster-cluster 
mergers play a prominent role in shaping the properties of cluster cores, in particular 
causing departures from hydrostatic equilibrium, and possibly disturbing cool cores. 
Our results do not support recent claims that large Einstein radius clusters present a 
challenge to the CDM paradigm. 

Key words: Gravitational lensing - Galaxies: clusters: general - Galaxies: clusters: in- 
dividual (A521, A611, A773, A868, A1413, A1835, A2204, RXJ1720, RXJ2129, Z2701) 



1 INTRODUCTION 

The evolution of galaxy clusters with cosmic time is an im- 
portant cosmological probe, as it traces the gravitational 



growth of dark matter on large scales. In particular, the 
cluster mass function can be directly tested against cosmo- 
logical models, as it is related to fundamental parameters 
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such as the matter den sity S2 m , the normali s ation of the 
powe r spectrum as (e.e. ISchuecker et al.ll2003l : ISmith et al.l 
2003), and the dark energy equation of state parameter w. 

Commonly used probes of the mas s of clusters at z ~ 
0.1-0.5 include the A"-band luminosi ty ihin et al.ll200rj). X - 
ray luminosity and temperat ure (e.g. IVikhlinin et al.ll2006Tl. 
Sunayev-Zeldovich effect (e.g. lNagaill2006l;lBonamente et al.l 
l2008h . cluster kinematics (e.g^lBlindert et al.l 12004), and 
gravitational lensing (e. g. iBardeau et al.l 120071 : ISmith et al.l 
l2005l : lOkabe et al.i r2009) . Systematic errors in these various 
mass probes can be calibrated by an inter-comparison of 
results from the various methods. The pursuit of the most 
robust calibrations based on low redshift cluster samples is 
essential to achieve reliable results from future studies at 
higher redshift, and to c ontrol systematic un certanties in 
cosmological experiments l|Albrecht et al.ll2006l ). 

Gravitational lensing plays a central role in this effort, 
as it does not rely on assumptions about the symmetry or 
equilibrium properties of the cluster mass distribution, al- 
though parametrized lensed model implicitly assume, for in- 
stance, elliptical symmetry. Indeed, the highest quality mea- 
surements of the dark matter mass and substructure in low 
redshift clusters have been obtained through the combina- 
tion of strong and weak gravitational lensing. The identifi- 
cation of background galaxies forming strong lensing arcs in 
cluster cores is a direct measurement of the enclosed mass 
within the Einstein radius 9e, providing an accurate nor- 
malisation of cluster mass m odels both in the c ore, and ex- 
tending to larger radii (e.g. iKneib et al.ll2003l) . The main 
caveats on lensing studies are that the lensing signal is sen- 
sitive to the mass distribution projected along the line of 
sight through the cluster, that the commonly used paramet- 
ric models must by design assume a parametric form of the 
mass distribution, and that gravitational lensing may prove 
to be impractical in samples of more distant ( z > 1) cluster 
samples that are starti ng to be discovered (e.g. iRosati et al.l 
ll999l ; lKurk et al.ll2009l) . 

Early joint lcnsing/X-ray cluster studies identified large 
discrepancies betwe en total cluster masses obtained from the 
two methods (e.g. iMiralda-Escude fc Babull 119951 ). These 
discrepancies were subsequently eliminated, albeit within 
quite large uncertainties, in cool core clusters prov ided that 
a two phase gas model was used (|Allenl 1 19981 ; see also 
ISmail et al.l 1 1997T ) . At a similar time, HST data started to 
became available for some clusters, allowing more precise 
strong lensing models t o be constructed (e.g. IKneib et al.l 
1996; iTvson et al.|[l998T) . More recent strong lensing stud- 
ies of cluster mass distributions have generally continued 
to concentrate on detailed studies of spectacular individual 
strong lensin g clusters based on deep multi-filter HST obser- 
vations (e.g. Broadhurst et al.1120051 : iLimousin et al.ll2007bl : 
iRichard et alT 2009) . The interpretation of results from such 
studies, e.g. size of Einstein radii and shape of density pro- 
file, in the context of the general cluster population is in- 
evitably proble matic. 

In contrast. ISmith et al.l l|2005l . hereafter Sm05) studied 
10 X-ray selected clusters (I/jc,o.2-2.4koV > 4.1 x 10 44 ergs _1 ) 
at z ~ 0.2 as a first step towards building large statistical 
samples of strong lensing clusters. Sm05 combined strong 
and weak lensing constraints from moderate depth HST ob- 
servations with X-ray spectro-imaging from Chandra, and 
near-infrared photometry of cluster galaxies from UKIRT. 



The main results included the first mass-observable scaling 
relations to employ lensing-based mass measurements for a 
well-defined sample; the only previous exampl e employed a 
comp ilation of 6 clusters from the literature |Hiorth et al.l 
1998). The main limiting factor on Sm05's results were the 
small sample size, with just 5 of the 10 clusters containing 
spectroscopically confirmed strongly lensed galaxies. 

The Local Cluster Substructure Survey (LoCuSS, PI: 
G. P. Smith) is an all-sky systematic survey of 165 X- 
ray luminous clusters at 0.15 < z < 0.30 selected from 
the ROSAT All-sky Survey catalogs (|Ebeling et all Il998l . 
l200d : iBohringer et a.1.1 [2004] ) . In addition to seeking to im- 
prove the statistical precisi on by enlarging the s ample size 
of studies such as Sm05 and lBardeau et aL (|2007l ). LoCuSS 
aims to incorporate new constraints on cluster baryons, most 
notably from observat ions of the Sunyaev-Zeldovich Effect 
l|Carlstrom et all 120021) to deliver a definitive local calibra- 
tion of mass-observable scaling relations that will be useful 
for cluster cosmology, interpretation of high-redshift cluster 
samples, and probing the physics of gas heating and cooling 
in merging clusters at low redshift. More generally, LoCuSS 
aims to constrain the scatter in the baryonic properties of 
the local cluster population and to correlate the scatter with 
the recent hierarchical assembly history of the clusters, the 
latter being constrain ed by the lensing-based mass maps 
(|Smith fc Tavlorll2008l ). A wide variety of studies are there- 
fore underway, including wid e-field weak-lensing analysis 
with Suprime-Cam on Subaru (lOkabe et al. 200S), th e first 
lensing-based mass-Ysz relation (|Marrone et al. 2009), and 
multiwavelength studies supported by space-based (HST, 
Herschel, Chandra, GALEX, Spitzer) as well as ground- 
based (Keck, VL T Gemini, MMT CTIO, KPNO, Pal o- 
mar) facilities fe.g.lZhang et al.ll2008l:lHaines et al.ll2009al lbl: 
ISanderson et al.ll2009bl ; ISmith et al.ll2009bl) . 

This paper presents a four-fold increase in the number 
of clusters with spectroscopically confirmed strong lensing 
clusters on Sm05. Such constraints remain critical within 
the overall context of the studies outlined above because 
strong-lensing offers precise constraints on the mass dis- 
tribution in cluster cores (R ~ lOOkpc) which for exam- 
ple provide an invaluable constraint on small scales when 
trying to measure the shape of cluster density profiles in 
conjunction with wide-field weak-lensing data from Subaru. 
We compare the details of the central mass distributions 
(integrated mass and substructure fraction) obtained from 
the strong-lensing constraints with near-infrared and X-ray 
probes of the cluster baryons. We describe the data in ij2] 
and present the strong-lensing analysis and models in Sj3] 
The main results on the mass and structure of the cluster 
cores are presented in Sj4] and summarized in Sj5] Through- 
out the paper, we use magnitudes quoted in the AB system, 
and a standard A-CDM model with Q m = 0.3, Qa = 0.7, 
and Ho = 70kms~ 1 Mpc~ 1 , whenever necessary. In this cos- 
mology, 1" is equivalent to 3.3kpc at z = 0.2. We adopt the 
definition e = 1 — b/a of the ellipticity, where a and b are the 
semi-major and semi-minor axis of the ellipse, respectively. 



2 OBSERVATIONS AND DATA REDUCTION 

We present here the cluster sample and the relevant datasets 
used in our study. The strong lensing analysis relies on two 
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Table 1. Optical and near-infrared imaging observations of the new strong-lensing LoCuSS clusters. 
Cluster ctBCG ^BCG 2 HST Near-infrared 











PID 


Camera/Filter 


Depth (Za) 


Camera 


J(Za) 


K s {Za) 


K* 

b 


A521>) 


73.528753 


-10.223605 


0.2475 


11312 


WFPC2/F606W 


27.0 


WIRC 


21.5 


21.0 


17.60 


A611 


120.236680 


36.056725 


0.2850 


9270 


ACS/F606W 


27.7 


WIRC 


21.0 


20.9 


17.85 


A 773 


139.472660 


51.727024 


0.2170 


8249 


WFPC2/F702W 


27.6 


WIRC 


21.7 


21.5 


17.36 


A 868 


146.359960 


-8.651994 


0.1535 


8203 


WFPC2/F606W 


27.5 


ISPI 


21.9 


21.6 


16.70 


Z270l( b ) 


148.204560 


51.885143 


0.2140 


9270 


ACS/F606W 


27.7 


WIRC 


21.5 


21.5 


17.33 


A 1413 


178.824510 


23.404451 


0.1427 


9292 


ACS/F775W 


27.3 


WIRC 


23.0 


22.3 


16.55 










9292 


ACS/F850LP 


26.7 










A 1835 


210.258860 


2.878532 


0.2528 


8249 


WFPC2/F702W 


27.7 


ISPI 


22.0 


21.3 


17.64 










10154 


ACS/F850LP 


27.3 










A 2204 


248.195540 


5.575825 


0.1524 


8301 


WFPC2/F606W 


26.4 


FLAM. 


22.1 


21.9 


16.68 


RXJ1720< C ) 


260.041860 


26.625627 


0.1640 


11312 


WFPC2/F606W 


26.9 


WIRC 


20.6 


20.7 


16.82 


RXJ2129M 


322.416510 


0.089227 


0.2350 


11312 


WFPC2/F606W 


26.7 


WIRC 


21.7 


21.7 


17.51 



M also known as RXC J0454.1-1014 00 also known as ZwCl 0949.6+5207 < c > RX J1720.1+2638 < d > RX J2129.6+0005 



Table 2. Extended sample of 10 previously published spectroscopically-confirmed strong-lensing LoCuSS clusters 



Cluster «BCG ^BCG z Strong-lensing reference 



A 68 


9.278626 


9.156722 


0.2546 


Sm05. Richard et al. (2007): Smith et al. ( 20021 


A 383 


42.014079 


-3.529040 


0.1883 


Sm05. Smith et al. (20011; Sand et al. (20041 


A 963 


154.264990 


39.047228 


0.2050 


Sm05. Ellis et al. (19911 


A 1201 


168.227080 


13.435946 


0.1688 


Edge et al. (20031 


A 1689 


197.872950 


-1.341005 


0.1832 


Limousin et al. (2007b1 


A 1703 


198.771971 


51.817494 


0.2800 


Richard et al. (20091 


A 2218 


248.954604 


66.212242 


0.1710 


Sm05. Kneib et al. (19961. Eh'asdottir et al. (20071 


A 2219 


250.082380 


46.711561 


0.2281 


Sm05 


A 2390 


328.403290 


17.695740 


0.2329 


Julio, PhD thesis 


A 2667 


357.914125 


-26.084375 


0.2264 


Covone et al. (20061 



main ingredients: (a) the identification of multiply-imaged 
systems and (b) measurements of their spectroscopic red- 
shifts. The former relies on high angular resolution HST 
imaging, while the latter is carried out using sensitive optical 
multi-object spectrographs on large aperture ground-based 
telescopes, in this case LRIS on the Keck-I 10-m telescope. 
We also use near-infrared photometry to select likely clus- 
ter galaxies based on the prominent red sequence of cluster 
galaxies in the J — K/K color-magnitude diagrams, that is 
largely insensitive to spectral type at z ~ 0.2. These photo- 
metric catalogs are also used in the paper during the con- 
struction of the lens models ( £13, 2 [I . Finally, we summarize 
the analysis of X-ray data available for the majority of these 
clusters, and which we latter use as a comparison with the 
strong-lensing results (see Sj4]). 



2.1 Cluster sample 

The sample in this article comprises the 10 clusters that 
we have spectroscopically confirmed as strong-lenses during 
our Keck observing campaign 2004-2008 (Table [1]), plus a 
further 10 clusters within the LoCuSS sample that have been 
previously published (Table . To qualify for observations 
with Keck, a cluster must have been previously imaged with 
HST (79 clusters), and lie at 8 > -25°. We note that some of 
the clusters with new strong-lensing models had previously 
been identified as candidate strong-lensing clusters: 



• A 773 and A 1835 were included in Sm05's sample how- 
ever in the absence of spectroscopic confirmation, only weak- 
lensing constraints were u sed to construct the mass m odels. 

• A 521 was studied by iMaurogordato et al.l (|2000h using 
ground-based images. A giant arc feature was observed but 
spectroscopic follow-up attempts were unable to confirm its 
strong lensing origin. 

• A 773, A 8 6 8, A 1835, A 2204 and Z2701 were part of 
the lSand et all (|2005h sample, who performed a systematic 
search for radial and tangential arc features in archival HST 
images of massive clusters. They identified radial arcs in 
A 773, A 1835, and A 2204. 

2.2 Imaging data 

2.2.1 HST imaging data 

High resolution imaging data taken with the ACS or 
WFPC2 instrument on HST are available for each selected 
cluster in one or two bands, either through our dedicated 
LoCuSS program (GO-DD 11312, PI: G. P. Smith), or from 
the archive. These data are summarized in Table [1] In every 
case, one of the observed filters is located between the V 
and the I band, the most common being the F606W filter. 

Each HST observation consists of several independent 
frames arranged in a dither pattern. These images were com- 
bined using the multidrizzle software with a final plate 
scale of 0.05" (ACS) or 0.1" (WFPC2). In the case of A 868, 
the observations comprise a mosaic of 6 WFPC2 pointings 
taken at different epochs. This mosaic was combined using 
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the relative astrometric shifts measured from bright objects 
in common between two adjacent pointings. 

A careful visual inspection was performed on the re- 
duced HST images by two of the authors (JR and GPS) 
to search for strong lensing features. A catalog of candidate 
multiple images was compiled by selecting highly sheared 
arc-like features, as well as objects displaying a very sim- 
ilar morphology and the characteristic symmetric effect of 
strong lensing (e.g., Fig. [2]). The majority of these multiple 
images does not appear like arcs, but have a clumpy mor- 
phology, with multiple knots of star-formation. These lensed 
sources are not d etected by automatic arc fin ding techniques 
(e.g. Arcf inder. ISeidel fc Bartelmanr] [2007). and therefore 
justify the use of a visual inspection of the images. 

2.2.2 Near-infrared imaging data 

J- and its-band data were obtained between March 2003 
and April 2007 on the following near-infrared instruments: 
WIRC on the Palomar-200in telescope, ISPI on the CTIO 
Blanco 4-m telescope, and FLAMINGOS on the Kitt Peak 
(KPNO) 4-m telescope. The properties of the data are sum- 
marized in Table [1] Observations in each filter comprised 
multiple frames with individual exposure times in the range 
30 — 120sec. These images were combined using standard 
IRAF reduction techniques, full details will be provided in 
Smith et al. (in preparation), and flux-calibrated using the 
2MASS source catalog. A dithering box of 80" width was 
used in order to improve the photometry of the extended 
envelope of the BCGs. The field of view of the final dithered 
images is larger than 8x8 arcmin 2 , centred on the Bright- 
est Cluster Galaxy (BCG). This entirely covers the central 
R < 600 kpc of each clus ter (as measured from t he BCG). 
We used the SExtractor (jBertin fc Arnouts 1996) software 
to derive a photometric catalog of MAGAUT0 magnitudes in 
the Ks band, as well as the J — Ks colours inside 4.0- 
"diameter apertures. 

2.2.3 Photometric catalogs of cluster galaxies 

In order to include galaxy-scale components in the strong- 
lens models ( i]3.2|l . it is important to perform a systematic 
selection of cluster members with their geometrical parame- 
ters (centroid a c , 5 C , ellipticity e c and position angle 9 C ) and 
total photometry in the Ks band. Indeed, Ks magnitudes 
are a more accurate proxy for stellar mass when making the 
assumption that light traces mass. 

Cluster galaxies were selected from the near-infrared 
catalog using the red seque nce technique on a J — Ks/K s 
colour-magnitude diagram llVisvanathan fc Sandagdll977h . 
We adopted a liming magnitude of Ks = 20 (MAGAUTO) and 
a colour width of 0.3 magnitudes for the sequence (above 
3 times the photometric uncertainties). The corresponding 
diagrams and colour selection regions are presented in Fig. 
[T] A final catalog of cluster members is obtained by corre- 
lating the HST and near-infrared photometric catalogs: in 
the central region, the geometrical parameters (a c , 8 C , e c , 
9 C ) of cluster members are replaced by their more accurate 
measurements from the HST image. 

As a sanity check, we plot the red sequence parameters 
used in our photometric selection (J — Ks colour at fixed 



its-band magnitude, and linear slope kjk s of the colour- 
magnitude diagrams) as a function of the cluster redshift 
(Fig. [U last 2 panels). The correlation between J — Ks 
at Ks = 18 follows the expected co lour of elliptica l galax - 
ies (using the empirical model from IColeman et al.l (|l980l ). 
while the slope kjk s is scattered around an average value 
kjk s = —0.028, compatibl e with similar measurements ob- 
tained bv lStott et all (2009) in this range of cluster redshifts. 

We estimated Ks, the Ks ba nd magnitude of an L* 
galaxy, following the results from iLin et all (120061) . They 
show that the evolution of K s with cluster redshift is best 
reproduced by a single b urst model at z/ orm = 1.5 from 
iBruzual fc Charlotl (|2003l) , assuming a Salpeter initial mass 
function and solar metallicity. Their best fit K s (z) values are 
reported in the final column of Table[T] Our current selection 
down to Ks = 20 enables us to select cluster members down 
to ~ 0.1L*. Since the fraction of fainter cluster members not 
included in our selection contribute to less than 1% of the 
total it-band luminosity, this selection does not affect our 
results on the it-band properties in the entire sample. 

2.3 Spectroscopic data 

2.3.1 Observing strategy 

We us ed the Low Re solution Imager and Spectrograph 
(LRIS, lOke et al.|[l995h on the Keck-I telescope to perform 
long-slit and multi-slit observations of the clusters. The spec- 
troscopic data used in the current paper is the outcome of 6 
different observing runs between 2004 and 2008, which are 
summarized in Table [3] 

We designed multi-slit masks containing ~ 30 slits of 
1.0" width to include as many of the candidate multiple 
systems as possible, with a few tilted slits following the ge- 
ometry of long arcs. In the case of fainter or less reliable iden- 
tifications, several images of the same system were observed 
in separate slits. Additional slits in the mask included faint 
galaxies with unknown redshifts, bright infrared sources se- 
lected at 24/im, or cluster members. 

The dichroic/grism/grating combinations listed in Ta- 
ble [3] to ensure a high throughput over the wavelength range 
3500 to 9500 A and good spectral resolution (Table[3| in the 
red to resolve OH skylines as well as [Oil] emission line dou- 
blets in the redshift range 0.85 — 1.50. The total exposure 
time varied between 1 and 3 hours per mask, depending on 
the seeing conditions and the detection of strong emission 
lines when looking over a first reduction of the data at the 
telescope. 

2.3.2 Data reduction 

The data we re reduced using the Python version of the 
iKelsonl l)2003h reduction scripts, which offer the advantage 
of processing the images in their distorted framework. This 
helps to reduce noise correlations, in particular for the case 
of tilted slits. We performed standard reduction steps for 
bias removal, flat-field correction, wavelength calibration, 
sky subtraction, and cosmic-ray rejection, and used observa- 
tions of spectroscopic standard stars obtained on the same 
nights to derive the flux calibration. 

The extracted spectra of objects used in the present 
paper are shown in Fig. IA1I (appendix) . These spectra show 
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Figure 1. J — Kg vs Kg colour-magnitude diagrams for each observed cluster, showing the red sequence selection of cluster members 
(red dots). The photometric selection used is shown as delimited by the solid box (see text for details). The last two panels show the red 
sequence parameters (intercept at Kg = 18 and slope kjk s ) as a function of cluster redshift. The solid line gives the expected colours 
of an elliptical galaxy, and the dashed line marks the average slope kjk s f° r the sample. 

A 1413 z = 0.1427 A 2204 z = 0.1524 A868 2 = 0.1535 




RXJ2129 z = 0.2350 A 521 z = 0.2475 A 1835 z = 0.2528 




" 1 1 1 1 1 1 1 1 1 1 0.1 0.15 0.2 0.25 0.1 0.15 0.2 0.25 

12 14- 16 18 20 



Table 3. Summary of the Keck/LRIS spectroscopic data. 



Cluster 


Date 


Config 




Grism 


Grating[A c ] 


Dichroic 


(ks) 


Seeing 


Standard star 


A 521 


Feb 2007 


MOS 


18 


400/3400 


900/5500 


[6320] 


560 


3.8 


1.1" 


Feige 92 










400/3400 


831/8200 


[8100] 


560 


3.6 


1.1" 


Feige 92 


A611 


Nov 2006 


LS 


1 


600/4000 


400/8500 


[7970] 


560 


4.5 


0.9" 


Feige 110 




Jan 2007 


LS 


1 


600/4000 


400/8500 


[7970] 


560 


1.8 


1.3" 


Feige 34 


A 773 


Feb 2007 


MOS 


26 


400/3400 


900/5500 


[6400] 


560 


3.6 


1.0" 


Feige 92 










400/3400 


831/8200 


[8070] 


560 


3.6 


1.0" 


Feige 92 


A 868 


May 2008 


MOS 


36 


300/5000 


600/7500 


[8100] 


680 


5.4 


0.7" 


BD+33-2642 


Z2701 


May 2008 


MOS 


22 


300/5000 


600/7500 


[8100] 


680 


3.6 


1.0" 


BD+33-2642 


A 1413 


Feb 2004 


LS 


1 


600/4000 


400/8500 


[7640] 


560 


5.4 


0.9" 




A 1835 


March 2005 


MOS 


36 


400/3400 


600/7500 


[6850] 


560 


7.2 


0.7" 


G138-31 


A 2204 


May 2008 


MOS 


37 


300/5000 


600/7500 


[8100] 


680 


3.6 


0.7" 


BD+33-2642 


RXJ1720 


May 2008 


MOS 


21 


300/5000 


600/7500 


[8100] 


680 


5.4 


0.9" 


BD+33-2642 










400/3400 


600/7500 


[8100] 


680 


5.4 


0.9" 


BD+33-2642 


RXJ2129 


May 2008 


MOS 


22 


300/5000 


600/7500 


[8100] 


680 


5.4 


1.1" 


BD+33-2642 
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Table 4. Properties of the multiple imaged systems. Redshift values quoted with brackets are predictions from the lensing model. The 
last two columns give the observed HST magnitudes and linear magnification factor predicted by the lens model. 



Cluster / Multipl 


e a 


S 


z 


Reference 


HST 


fx (mags) 


Lt (linear) 


A 521 










(F606W) 






1.1° 


73.526932 


-10.223421 


1.043±0.01 


[Oil], CI I, Fell 


21.49 ±0.04 


1.76±0.10 


5.1±0.5 


1.2 


73.527544 


-10.222524 


1.043±0.01 


CI I, Fell 


21.32 ±0.04 


1.71±0.10 


4.8±0.4 


1.3 


73.529391 


-10.221431 


1.043±0.01 


[Oil] 


21.57 ±0.20 


0.98±0.10 


2.5±0.2 


A611 










(F606W) 






1.1° 


120.23219 


36.061593 






24.58 ±0.03 


1.38±0.11 


3.6±0.4 


1.2 


120.24171 


36.055210 


2.06 ±0.02 


cm, CIV 


23.87 ±0.05 


2.52±0.30 


10.2±2.8 


1.3 


120.24103 


36.058273 


2.06 ±0.02 


cm, civ 


22.80 ±0.10 


1.95±0.18 


6.0±1.0 


1.4 


120.23555 


36.054247 






24.98 ±0.04 


0.79±0.06 


2.1±0.1 


1.5 


120.23589 


36.054887 






25.46 ±0.05 


1.0±0.07 


2.5±0.2 


2.1 


120.23716 


36.061183 


0.908 ±0.005 


[Oil] 


23.26 ±0.02 


2.47±0.29 


9.7±2.6 


2.2 


120.24046 


36.059825 


0.908 ±0.005 


[Oil] 


23.26 ±0.09 


1.10±0.08 


2.8±0.2 


2.3 


120.24206 


36.057539 


0.908 ± 0.005 


[Oil] 


23.17 ±0.14 


2.48±0.29 


9.8±2.6 


3.1 


120.23554 


36.060897 


[2.5t^] 




24.42 ±0.03 


4.08±1.27 


42±15 


3.2 


120.23739 


36.060653 




25.05 ±0.06 


2.04±0.19 


6.5±1.2 


3.3 


120.24311 


36.053638 






25.94 ±0.08 


1.38±0.11 


3.6±0.4 


3.4 


1Z0.ZO4U t 


36.0558U/ 






25.63 ± o.u / 


i.3y±u.ii 


3.6±0.4 


4.1 


120.24182 


36.056242 


2.59 ±0.01 


Ly a 


27.14 ±0.19 


4.35±1.62 


55±17 


4.2 


120.23193 


36.062088 






27.67 ±0.19 


1.26±0.09 


3.2±0.2 


A 773 










(F702W) 






1.1 


139.48925 


51.729656 


2.300±0.005 


Ly a , Sill, SUV, CIV 


22.95 ±0.02 


2.69±0.35 


11.9±3.8 


1.2 


139.48920 


51.724870 


2.300±0.005 




23.03 ±0.02 


3.0±0.5 


16.0±6.9 


1.3 


139.48351 


51.741593 






23.47 ±0.02 


1.57±0.13 


4.2±0.5 


2.1 


139.46920 


51.732350 


3.84±0.01 


Ly a 


23.85 ±0.06 


2.42±0.27 


9.3±2.3 


2.2 


139.48762 


51.716864 






23.58 ±0.06 


1.18±0.09 


3.0±0.3 


3.1 


139.49535 


51.728050 


1.010±0.005 


[Oil] 


23.09 ±0.05 


2.0±0.2 


6.1±1.0 


3.2 


139.49501 


51.728985 


1.010±0.005 


[Oil] 


21.54 ±0.05 


2.0±0.2 


6.1±1.0 


A 868 










(F606W) 






1.1 


146.36459 


-8.6481919 


0.551±0.002 


[OII],[OIII] 


22.60 ±0.03< 6 ) 


3.29±0.61 


36±7.8 


1.2 


146.36474 


-8.6490159 


0.551±0.002 


[OII],[OIII] 


22.60 ±0.03< 6 ) 


2.28±0.24 


19±4.7 


1.3 


146.36554 


-8.6502853 


0.551±0.002 




23.50 ±0.05 


2.64±0.34 


23. ±4.1 


Z2701 










(F606W) 






1.1 


148.20940 


51.885620 


1.163 


[Oil] 


23.73 ±0.03 


5.32±0.98 


20.7±6.1 


1.2 


148.20708 


51.883841 






24.22 ±0.03 


3.90±0.88 


8.2±1.8 


1.3 


148.20087 


51.882325 


1.163 


[Oil] 


24.34 ±0.03 


2.62±0.33 


11.4±3.7 


2.1 


148.21056 


51.884031 


P-5±8:§] 




25.05 ±0.08 


3.73±0.92 


130±60 


2.2 


148.20988 


51.883375 




24.99 ±0.08 


3.90±0.88 


36. ±9 


2.3 


148.20083 


51.880880 






25.52 ±0.05 


2.62±0.33 


11.2±3.4 


A 1413 










(F775W) 






1.1 


178.82496 


23.412185 


2.726±0.003 


Ly a , Sill 


23.16 ±0.06 


3.73±0.92 


31. ±6. 


1.2 


178.82328 


23.412089 


2.726±0.003 


Ly a , Sill 


22.79 ±0.08 


5.00±2.95 


100±27 


2.1 


178.82855 


23.399205 


2.030±0.004 


Ly a , Sill 


22.90 ±0.16 


2.37±0.26 


8.9±2.1 


2.2 


178.82591 


23.398715 


2.030±0.004 


Ly a , Sill 


22.99 ±0.23 


2.94±0.44 


15±6 


2.3 


178.81985 


23.399441 






23.82 ±0.09 


2.50±0.30 


10±2.8 


3.1 


178.82933 


23.406898 


1.20±0.01 


[Oil] 


23.03 ±0.06 


1.28±0.10 


3.3±0.3 


3.2 


178.82391 


23.407396 


1.20±0.01 


[Oil] 


23.64 ±0.21 


1.54±0.12 


4.1±0.5 


3.3 


178.81981 


23.405955 


1.20±0.01 


[Oil] 


22.84 ±0.13 


1.33±0.10 


3.4±0.3 


4.1 


178.82940 


23.409414 


[2.9±°;1] 




25.47 ±0.05 


1.30±0.10 


3.3±0.3 


4.2 


178.82263 


23.409657 




24.50 ±0.03 


2.07±0.20 


6.7±1.2 


4.3 


178.82007 


23.408814 






24.57 ±0.03 


1.95±0.18 


6.0±1.0 



location of the brightest knot in the image. O Photometry is the sum of 1.1 and 1.2 images 



Lyman-a (in absorption or emission) and/or additional ul- 
traviolet absorption lines, or the resolved [Oil] doublet. For 
some faint multiple systems observed in different slits, red- 
shift measurements were derived after stacking the relevant 
exposures. The average redshift value is obtained from the 



peaks of the main spectral features identified, while the cor- 
responding error is taken from the spectral dispersion. Ad- 
ditional uncertainties generated by the accuracy of the rela- 
tive and absolute wavelength calibrations, about 1.1 and 1.5 
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Cluster / Multiple 


Q 


5 


z 


A 1835 








1.1 
1.2 


210.26573 
21 26403 


2.8706065 
2 86Q3007 


[2.5l° . 2 2 ] 


1.3 


210.24740 


2.8698812 




2 1 

9 9 
z .z 


9i n 9fifii 4 

ZlU.ZOOOo 


9 fi71 /1 1 91 
Z .0 / 141zl 


11 3+ - 1 ! 


2.3 


210.24824 


2.8717531 




3.1 


210.26326 


2.8849801 




3.2 


210.26278 


2.8852175 


3.3 


210.24374 


2.8775381 




A 1 
1.1 


ZlU.zOooZ 


Z.oo4u / Oo 


[2-3_ J 


1 2 


91 fl 9^0.73 


9 SP^R^I Q 
z.oouooiy 


4.3 


210.24435 


2.8781489 




5.1 


210.25984 


2.8824075 


2.6 


5.2 


ni n c\ a A o A 

210.24484 


2.8721518 


2.6 


D . O 


91 fl 9^091 






6.1 
6.2 


210.26165 
210.26131 


2.8775836 
2.8777493 


[3-6I 14 ] 


7.1 


210.25395 


2.8731863 


2.070±0.004 


7.2 


210.25421 


2.8803387 




7.3 


210.27108 


2.8801023 




A 2204 








1.1 


248.19665 


5.5781303 


1.06±0.01 


1.2 


248.19457 


5.5689343 


1.06±0.01 


1.3 


248.19603 


5.5779196 




RXJ1720 








1.1 


260.04314 


26.624351 


2.136±0.005 


1.2 


260.04247 


26.624155 


2.136±0.005 


1.3 


260.03779 


26.626930 




1.4 


260.04285 


26.626111 





Reference 



HST 


/i (mag 




fi (linear 


(F850LP) 








O/l OO I n 10 


z.loztu 


99 
ZZ 


7 94-1 c; 


oo ciA _i_ n on 


4.oo±l 


OU 


er /I _l_ 1 c 


25.12 ± 0.25 


0.83±0 


06 


2.1±0.1 


Zo.Zo ± 0.2o 


1.60±0 


13 


4.4±0.5 


22.87 ± 0.29 


0.42±0 


1 


1.5±0.05 


oc AO 4- n 


0.85±0 


07 


2.2±0.1 


25.63 ± 0.07 


2.31±0 


25 


8.4±1.9 


ZO.ZO IE U.UD 


Z.OOIEU 


1 


ll.Dxj.D 


25.87 ± 0.07 


1.87±0 


17 


5.6±0.9 


24.70 ± 0.04 


1.96±0 


18 


6.1±1.0 


25.41 ± 0.07 


3.15±0 


5 1 


1 A 1 O f\ 

17.4±8.0 


25.08 ± 0.13 


0.83±0 


06 


2.1±0.1 


23.44 ±0.03 


0.67±0 


05 


1.8±0.1 


23.90 ± 0.02 


0.69±0 


06 


1.9±0.1 


26.04 ± 0.18 








25.77 ±0.09 


0.99±0 


07 


2.5±0.2 


25.43 ±0.06 


0.80±0 


06 


2.1±0.2 


22.28 ±0.04 


4.86±2 


59 


88±30 


23.96 ± 0.05 


4.00±1 


18 


40±21 


25.02 ± 0.05 


0.86±0 


07 


2.2±0.1 


(F606W) 








23.60 ±0.12 


1.25±0 


09 


3.2±0.3 


22.61 ± 0.03 


0.85±0 


06 


2.2±0.1 


23.70 ±0.15 


0.53±0 


05 


1.6±0.08 


(F606W) 








22.74 ± 0.05 (c) 


3.16 ±0.54 


18.4 ±9. 


22.74 ± 0.05 (c) 


1.54±0.12 


4.1±0.5 


24.85 ± 0.08 


0.91±0.07 


2.3±0.2 


24.41 ± 0.07 


1.33±0.10 


3.4±0.3 


(F606W) 








23.41 ± 0.03 


4.30±1.56 


52±17 


23.05 ± 0.02 


2.33±0.25 


8.6±2.0 


24.74 ± 0.06 


2.23±0.23 


7.8±1.6 



RXJ2129 



Ly a (a) 
Ly a (a) 



1.1 


322.42040 


0.088305 


1.2 


322.42018 


0.089722 


1.3 


322.41798 


0.093250 



1.965±0.005 SUV, Sill, CIV 
1.965±0.005 SUV Sill, CIV 



( a ) Dannerbauer et al., in preparation W See also lWilman et al.l d2006t) ( c ) Photometry is the sum of 1.1 and 1.2 images 



A, respectively, were quadratically added to yield the final 
redshift errors. The redshifts are listed in Table [4] 



2.4 X-ray observations 

X-ray properties of the LoCuSS clusters have been deter- 
mined for 18/20 clusters from Chandra archival observa- 
tions, which we r e redu ced and analyzed as described in 
ISanderson et al.l (l2009ah follow ing the methods developed 
bv lSanderson fc Ponmanl |2009l ). Annular gas temperature 
and density profiles were derive d and fitted with the phe - 
nomenological cluster model of lAscasibar fc Dieed {2008), 
to determine th e total clu s ter m ass profile. These models 
are based on the lHernauistl (|l990t) profile, which has the ad- 
vantage of yielding convergent three- and two-dimensional 
mass measurements. The latter is particularly helpful when 
comparing lensing and X-ray based measurements of the 
projected mass within R < 250kpc (see i j2.4l) . 

The X-ray models were also used to measure the gra- 
die nt of the logarithmi c gas density profile at 0.04rsoo 
(a; IVikhlinin et al.l l2007t> . which quantifies the strength of 
coolin g in the core l|Vikhlinin et al] l2007t ISanderson et al.l 
l2009al) . Core-excluded mean temperatures for the clusters 



were also measured in the r ange 0.15 — 0.2r5 o o, acco rding 
to the procedure outlined in I Sanderson et all 1 2009bl) . We 
classify clusters as hosting a cool core if the slope of the gas 
density profile is a < —0.85. This matches the range of a 
displayed by clu s ters th at contain an Ha emitting BCG in 
Sanderson et al ] i|2009al ). Under this definition, 8/18 clus- 
ters as cool core clusters, corresponding to 44 ± 14 % where 
the er ror bar is a binomial uncertainty following iGehrelsl 
(1 19861 ). We also classify clusters as "disturbed" if the offset 
between the X-ray centroid and the BCG is > O.Olrsoo, this 
identifies all except 3 of the clusters c ontaining an Ha emit- 
ting BCG in ISanderson et al.l (|2009al . see their Fig. 5). This 
yields a fraction of 67~^_\^% of the sample that are classified 
as undisturbed. All cool core clusters are also classified as 
undisturbed, however 4 non-cool core clusters are also classi- 
fied as undisturbed. The resulting X-ray properties are given 
in Table El 



3 STRONG-LENSING ANALYSIS 

We now describe the multiply-imaged galaxies used to con- 
strain mass models of the 10 new strong-lensing clusters 
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Table 5. Summary of Chandra X-ray properties for the Keck sample (1st half) and the extended sample (2nd half) 



Cluster 


L x 


ID 6 


Temperature c 
(keV) 


^500 
(kpc) 


M X (R < 250kpc) 
(10 13 Mq) 


a d 


Cool 
core? 


Offset e 
(kpc) 


Disturbed? 


A521 


9.45 


901 


7.08 + ^g3 


921 ± 26 


8.5 ±0.5 


-0.10 ±0.29 


No 


37 


Yes 


A611 


8.05 


3194 


7.m+lf 2 


1342 ± 140 


20.6 ± 1.9 


-0.70 ±0.04 


No 


1 


No 


A 773 


7.74 


5006 


7.50+? : T 


1358 ± 45 


18.0 ± 1.1 


-0.40 ±0.05 


No 


20 


Yes 


Z2701 


6.32 


3195 


_ 41 

5.08 n ' a 


1271 ± 300 


15.0 ±4.1 


-0.88 ±0.12 


Yes 


3 


No 


A 1413 


7.80 


5003 


6-90+n?? 

— U.OO 


1296 ± 28 


18.4 ±0.5 


-0.68 ±0.04 


No 


2 


No 


A 1835 


22.80 


6880 


9.82+H! 

— U.44 


1506 ± 57 


24.6 ± 1.0 


-1.17 ±0.05 


Yes 


5 


No 


A 2204 


12.57 


7940 


9.64+°'« 

— U.D4 


1420 ± 59 


23.8 ± 1.4 


-1.22 ±0.04 


Yes 


8 


No 


RXJ1720 


9.54 


4361 


7.96lg : g 


1349 ± 71 


19.2 ± 1.1 


-1.06 ±0.04 


Yes 


6 


No 


RXJ2129 


11.00 


552 


8.27tl:« 


1155 ±33 


16.1 ± 1.1 


-1.09 ±0.03 


Yes 


6 


No 


A 868 


3.46 


N/A 














A68 


8.81 


3250 


8 - 89 i3.37 


941 ± 56 


11.9 ± 1.8 


-0.25 ±0.02 


No 


52 


Yes 


A 383 


5.27 


2320 


5 01 + - 43 


1014 ±47 


12.1 ± 1.0 


-1.09 ±0.03 


Yes 


2 


No 


A 963 


6.16 


903 


ry yo + 0.52 

°"'°-0.77 


1210 ±39 


16.9 ±0.7 


-0.68 ±0.05 


No 


6 


No 


A 1201 


3.72 


4216 


5.56+°f 


1016 ±30 


10.7 ±0.6 


-0.65 ±0.13 


No 


11 


Yes 


A 1689 


16.27 


5004 


a OO + 0.68 
S -« b _0.88 


1451 ± 23 


25.7 ± 1.0 


-0.77 ±0.03 


No 


3 


No 


A 2218 


5.51 


1666 


7 ! 7+0.51 
'• i '-0.58 


1216 ±34 


16.0 ±0.6 


-0.33 ±0.03 


No 


11 


Yes 


A 2219 


12.07 


896 


n.52ir 9 


1786 ± 149 


26.9 ± 1.4 


-0.31 ±0.03 


No 


28 


Yes 


A 2390 


12.69 


4193 


Q 7Q+0-55 
»• ,o _0.56 


1437 ±68 


20.0 ± 1.0 


-0.94 ±0.02 


Yes 


2 


No 


A 2667 


15.78 


2214 


5 66+ 62 


1243 ± 115 


14.1 ± 1.6 


-0.89 ±0.05 


Yes 


3 


No 


A 1703 


8.66 


N/A 















a X-ray luminosity in the 0.1 - 2.4keV band from ROSAT llEbeling et al.l|l998tl200ol;lBohringer et al.ll2004l) . b Chandra observation 
ident if ier. c Measured between 0.15 and 0.2rnnn f see [ Sanderson et al. i l2009bh . d Logarithmic slope of the gas density profile at 0.04r5oo 
(see lVikhlinin et al.ll2007l: ISanderson fc Ponmanll2009f) . e X-ray centroid/BCG offset from lSanderson et all i2009al) . Errors are lcr. 



( H3.1|) . and the methods used to construct the models them- 
selves (SEH). 



3.1 Multiple images and spectroscopic redshifts 

The multiple-image constraints used for each cluster are de- 
scribed below, and summarized in Table 21 including their 
astrometry, photometry and available spectroscopic red- 
shifts. The locations of the multiple images are also marked 
in Fig. O The naming convention for multiple images in a 
given cluster is the following: (1.1, 1.2, 1.3, ... ) are individ- 
ual images in system 1, (2.1, 2.2, 2.3,...) in system 2, and so 
on. 



A 521 

The WFPC2 observations revealed a triply-imaged face-on 
spiral just 6" from the BCG, and containing many indi- 
vidual knots of star formation. The clear symmetry be- 
tween these knots (marked A to E in Fig. [2} confirms a 
strong-lensing system of 3 images,which we measure to be at 
z — 1.034 from [Oil] emission and UV absorption lines. We 
included the 3 brightest knots of each image as 9 individual 
constraints for th e mass modeling. Th e syste m was previ- 
ously identified bv lMaurogordato et al. ( 1996) using CFHT 
ground-based imaging. However the poor angular resolution 
of their data prevented them from identifying it as strongly- 
lensed, and subsequent sp ectroscopy was also inconclusive 
l|Maurogordato et al.|[200ol) . 



A 611 

We identify 4 systems, the first of which (1.1 to 1.5) clearly 
appears as 5 images of identical clumpy morphology, each 
of them composed of 5 individual knots (A to E) of star 
formation. We obtained a spectroscopic redshift for 1.2 and 
1.3, showing identical features of a z — 2.06 galaxy. The 
locations of the radial counter-images (1.4 and 1.5) are con- 
firmed by a preliminary mass model. The second system (2.1 
to 2.3) appears as a giant arc to the North-East of the BCG, 
for which we derive z = 0.908 from strong [Oil] emission. 
This giant arc is clearly affected by the presence of clus- 
ter substructure, with significant changes of curvature close 
to individual galaxies. The third system (3.1 to 3.4) was 
identified thanks to the symmetry between 3.1 and 3.2. The 
mass model predicts counter images that we identify as 3.3 
and 3.4. Finally, a very faint arc (4.1) serendipitously falling 
in one of the slits was confirmed at z = 2.54 with strong 
Lyman-a emission. The mass model confirmed it is located 
on the z = 2.5 critical line, and predicts the counter image 
that we identify as 4.2. 



A 773 

We identify 3 multiply imaged systems in this cluster. Sys- 
tem 1 is a triply-imaged system with a clear symmetry be- 
tween components 1.1 and 1.2. We derived a redshift z = 2.3 
for both images and identify the counter-image 1.3 with help 
from the mass model. The second system (2.1 and 2.2) is a 
radial image and counter-image at z = 3.84 from strong Ly a 
emission. Finally, we found a spectroscopic redshift z — 1.0 
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for the third system (3.1 and 3.2) curving around one of the 
cluster galaxies. 

A 868 

A giant arc containing three bright images (1.1 to 1.3) is 
clearly found around the second brightest galaxy to the 
north-east of the BCG. We derive its redshift at z — 0.551 
from the combination of [Oil] and [OIII] emissions. The 
curvature of this arc suggests it is strongly affected by a 
secondary mass component, as well as the main component 
centred on the BCG. 



RXJ1720 

This cluster is constrained by a giant arc (1.1 and 1.2) curv- 
ing around the BCG, with a spectroscopic redshift z — 2.136 
from the magnesium doublet. We identify two counter im- 
ages (1.3 and 1.4) with help from the mass model. 

RXJ2129 

Two compact sources with mirror symmetry (1.1 and 1.2) 
have been identified z = 1.965 from UV absorption lines, and 
the location of the counter image 1.3 was identified with the 
mass model. 



Z2701 

We identify a first system (1.1 to 1.3) as three compact im- 
ages ~ 10" from the BCG, confirmed to be at z — 1.163. 
A second system is formed by two faint symmetric arcs (2.1 
and 2.2), for which the mass model predicts a counter-image 
identified as 2.3. 



A 1413 

The two-filter ACS data enable 4 multiply-imaged galaxies 
to be identified from their colours and morphologies, three 
of them having spectroscopic redshifts between z = 1.17 
and 2 = 2.77, from [Oil] or Lyman-a in emission. System 
1 is a bright fold arc to the north, where we identify two 
bright knots as 1.1 and 1.2. The three other systems are 
triply-imaged with clear mirror symmetry. 



A 1835 

Several studies hav e previously identifi e d strongly-lensed 
arcs in this cluster l|Schmidt et al.ll200ll; ISand et al.l 120051 . 
Sm05) however to date spectroscopic redshifts have not been 
measured for any of these background galaxies, thereby lim- 
iting the precision of previous attempts to model the mass 
distribution in this cluster. We report redshift measurements 
for two systems: one of them is a triple system (labelled 5) 
discovered during a narrow-band search for Lyman-a emit- 
ters at z — 2.5 (Dannerbauer et al., in preparation), as two 
objects with similar morphologies and strong Ly-a emission. 
A radial counter image (5.3) is predicted close to the BCG. 
The second is a pair of merging images (system 7) measured 
at z = 2.07, for which we identify two counter images 7.2 
and 7.3. A 4th image 7.4 is predicted very close to the BCG 
but the light of the BCG prevents a reliable identification. 
Finally, we checked that the s ubmillimetre sources J1/J2 
discussed by ISmaife t al. (2005) at z — 2.56 are indeed pre- 
dicted by the mass model to be singly-imaged. 



A 2204 

The redshift of a radial arc (1.1 and 1.2) was measured 
to be z = 1.0 from strong [Oil] emission during an IFU 
observ ation of the central galaxy with the VIMOS instru- 
ment l|Wilman et al.l [20061 ) . We confirm this redshift inde- 
pendently, and find the same redshift for 1.3, which our mass 
model predicted to be the counter image. 



3.2 Cluster mass models 



We have used LenstooFl (|Kneiblll993l ; [jullo et al.ll2007t ) to re- 
construct the distribution of mass in the cluster cores, which 
is described as a superposition of analytic mass components 
that account for both cluster- and galaxy-scale mass, follow- 
ing Sm05. We use all the multiple-image systems described 
in the previous section and Table [4] as model constraints, 
using spectroscopic redshifts whenever available. The main 
differences between these models and those of Sm05 are 
that unknown multiple-image redshifts are free parameters 
in the models, and two different parameterisations are used 
for the smooth cluster-scale mass components ( H3.2.ip . and 
parameter space is explored using the Bayesian MCMC 
sampler that is now available within Lenstool ( H3.2.4p . In 
common with Sm05 and other studies of s t rong-lensing 
clust ers by our group, (e.g. IOi^ib_et_al. 1996): SjmthetjilJ 



200ll.l2002j; [Richard et al.ll2007l;lLimousin et al.ll2007bl . l2008l ; 
Richard et all 120091; ISmith et al.1 l2009al) ), the mass mod- 
elling is an iterative process, as we first produce a rudimen- 
tary mass model, based on the most reliable constraints, and 
then gradually incorporate additional constraints. 



3.2.1 Cluster-scale mass components 

The mass distribution of cluster cores within the strong- 
lensing region is generally well-described by a smooth ellip- 
tical profile centred near the BCG. As a starting point in our 
analysis, we try to reproduce the strong lensing constraints 
using a double pseudo-isothermal elliptical mass distribu- 
tion (dPIE, also know n as a smoothly truncated PIEMD, 
lElfasdottir et al.ll2007T ). The dPIE profile, which was also 
used in Sm05, is characterized by its central position, posi- 
tion angle and ellipticity, the central velocity dispersion 00 
and two characteristic radii: a core radius r corG and a cut 
radius r cut . In an (x, y) coordinate system oriented along 
the angle 9 of the ellipse, the surface mass density takes the 
form: 



T,{x,y) 



2G 



r C ut 



(rlo.e + p 2 ) 1 ^ (r c 2 ut + P 2 ) 1/2 



with p 2 = [(2 - e) (as - x c )/2] 2 + [(2 - e) (y - y e )/(2 - 2e)] 2 , 
and (x c ,y c ) is the centre of the mass distribution (see also 
Sm05). 



1 Publicly available at http: //www. oamp . f r/cosmology/lenstool 
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Figure 2. Multiple-image identifications and critical line at redshift 1.0 (A521), 0.91 and 2.1 (A611), 2.3 (A 773), 0.55 (A868), 1.2 
(Z2701) and 1.1 (A2204). Individual knots used as independent constraints are marked in A521 and A 611. 
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Continued. The critical lines are presented at redshift 2.07 (A 1835, system 7), 2.0 and 2.7 (A 1413), 2.1 (RXJ1720) and 2.0 (RXJ2129). 
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The mass models generally include a single dPIE com- 
ponent centred at the BCG location, but in some cases this 
description is unable to reproduce the observed multiple- 
image systems. When this is the case, we add a second 
or third cluster-scale mass component. A full discussion of 
our capacity to identify cluster-scale clumps from the strong 
lensing signal is provided in H3.4.3I 

In order to estimate the model dependence of our strong 
lensing results, we also independently fit the cluster-scale 
mass di stributions with the same number of Navarro-Frenk- 
White |Navarro et al.lll997l . hereafter NFW) profiles, char- 
acterized by a concentration C200 and a scale radius r s , de- 
scribed by usual density profile: 



p(r) 



PcSc 

(r/r a ) (1 + (r/r 3 )) 2 



(2) 



where p c is the critical density and S c is related to C200 
through the relation: 



200 



c 200 



3 ln(l + C200) - c 2 oo/(l + C200) ' 



(3) 



The generalisation of this spherical model into a projected 
elliptical mass distribution follows a similar relation as Eq.Q] 
The full analyt ical description of the pseudo-elliptical NFW 
is presented in lGolse fc Kneibl l|2002T ). 



3.2.2 Galaxy-scale mass component 

The presence of numerous cluster galaxies near the cluster 
centre, as identified by our near-infrared photometric selec- 
tion ( £12.2.31) affects the shape of the mass distribution lo- 
cally, and therefore the associated lensing signal. Following 
iKneib et al. ( 1996) and Sm05, we account for galaxy-scale 
mass in the models by adding individual dPIE mass compo- 
nents at the location of each galaxy in the cluster catalogs 
discussed in i]2.2.3l within 250kpc of the respective BCGs. 
The geometrical parameters of these components (x c , y c , e, 
8) are matched to the values measured by SExtractor. The 
parameters do, r COIC and r cu t, upon which the total mass of 
each galaxy-scale dark matter halo depend, are assumed to 
scale with their if -b and luminosity Lg, rel ative to an L* K 
galaxy, following the lFaber fc Jackson! 1 19761) relation and a 
constant mass-to-light ratio: 



'core — ?" corc (Lk/L k ) / , 

fcut = r* ut (L K /L* K ) 1/2 , 
cr = Oo {Lk/L* k ) 1/a 



(4) 



Similar scaling re l ations have been used in the past 
ilCovone et al.ll2006l : iLimousin et al.ll2007bl : lElfasdottir et al.l 



2007), showing that r* OTC is small and has little effect on the 
modell ing results. We therefo re adopt r* nrR = 0.15kpc, fol- 
lowinglLimousin et al.l l|2007bl ): see also lBrainerd fc Specianl 
i|2003h . 

The two remaining parameters (oq, r* ut ) are degen- 
erate, because the perturbation of the cluster lensing sig- 
nal caused by an individual galaxy depends on the total 
galaxy mass, given by a\ r cut . After performing various tests 
with Lenstool, we concluded that fixing one of these pa- 
rameters does not affect the ability of the models to re- 
produce the observed multiple-image systems, provided that 



galaxies located within 2" of a highly magnified arc are ex- 
cluded from the scaling relations (see 33.2.3[) . We therefore 
adopted a fixed r£ ut = 45k pc, which matche s recen t galaxy- 
galaxy weak-lensing results iNataraian et al.l (|2009l) . For our 
adopted value of r* ore = 0.15kpc the measured galaxy ve- 
locity dispersions (a) and dPIE velocity dispersions (cto) are 
in good agreement - 0.94 < a/ao < 1.04 - the range of val- 
ues arising from the ra nge of radii at which a is measured 
l)Elfasdottir et all l2007f > . We therefore emp loy a Gaussian 
prior o f q-q = fl58±27)kms~ 1 . making use of lBernardi et al.l 
(2003)'s observational results on a* . 



3.2.3 BCG and individual galaxy perturbers 

We model separately the brightest cluster galaxy (BCG) in 
each cluster using a similar dPIE profile as other cluster 
members, but optimizing the values of cto.bcg arm fcut.BCG 
independently. The same approach was used for 6 local 
galaxy-scale perturbers in A 611, A 773, A 1413, A 1835, 
identified in Fig. [2] 



3.2.4 Optimization 

The models are fitted to the multiple-image constraints us- 
ing the new Bayesian Markov chain Mo nte-Carlo (h e reafte r 
MCMC) sampler, described in detail in I Julio etaD (|2007h . 
This process uses the observational constraints (positions of 
the multiply imaged systems) to optimize the parameters 
describing the mass distribution by matching the location 
of each image of a given system in the source plane. The 
quality of the models can be estimated using the root mean 
square (RMS) deviation in the image-plane from the ob- 
served positions of the multiple-image positions predicted 
by the model, defined as: 



(xobsj.fc - xpred • fc ) 2 + (yobs - k - ypred 



(5) 



where (xobsj,^, yobs^ fc ) and (xpred^ fe , ypred^ fe ) are the 
observed and predicted location s of image .7 in system 
fc, respectively . Foll owing Sm05, ILimousin et all (|2007bD : 
IRichard et~alj < |2009h we have used 0.2" as the positional 
uncertainty of the multiple image identifications. 



3.3 Best fit parameters 

The best fit parameters of the mass models are listed in 
Table [|J] for both models that parameterize the cluster- 
scale mass components as dPIE and NFW profiles. The 
parameter uncertainties are based on the M CMC chains 
generated by Lenstool fsee lJullo et al.l 120071 . for details). 
We obtain <t; < 0.5" for 9/10 clusters, with 7/10 having 
<7j <C 0.2". These results are typical of strong lensing studies 
using a similar number of multiples images llRichard et al.l 
120071 : ILimousin et al.ll2007bl: IRichard et al.ll2009h . However 
in the case of A 1835, where we have the largest num- 
ber of multiple-image constraints (7 systems), we obtained 
Oi ~ 3.15", which is comparable w ith recent results on 
A 1689 (2.87". iLimousin et al.ll2007bh using 32 systems. We 
discuss the origin of this large value of Oi for A 1835 in H3.4.1I 
The best fit ellipticity e and position angle 6 of the 
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Table 6. Best-fit parameters of the mass models. For each mass component, we give the centre, ellipticity, orientation, core and cut 
radii, as well as central velocity dispersion of the dPIE profile. The following column gives the image plane RMS of this model, and the 
2 rightmost columns present the best-fit concentration and scale radius of the corresponding NFW profile, used to derive the systematic 
errors on the model. 
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cluster-scale dark matter halos are in good agreement be- therefore adopt the dPIE models as the fiducial models for 

tween the dPIE and NFW models, with < e < 0.7, and 6 the rest of the paper, and use the NFW models to quantify 

of the main cluster-scale halo agreeing with the orientation systematic errors in i)3.4l 

of the BCG. A few clusters have e > 0.5; |Golse fc KneiE For 8/10 clusters a single dPIE or NFW cluster-scale 

(|2002|) have shown that the pseudo-elliptical NFW profile dark matter halo is needed in addition to the cluster galax- 

presents a boxy /peanut shape at such high ellipticities. We les to reproduce accurately the observed multiply-imaged 
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thus dropping the assumption ot elliptical symmetry. The 
resulting mass map shows a mass extension to the North of 
the BCG and less extended to the South, compared to the 
original dPIE model. The quality of the fit is improved to 
Oi — 1.7" for the same multiply-imaged systems. The lack 
of spectroscopic redshifts in the Northern multiple images 
precludes strong conclusions about this asymmetry. Secur- 
ing spectroscopic redshifts for these multiple-image systems 
is therefore a necessary step to further improve this mass 
model. Overall, the absence of a secondary cluster-scale com- 
ponent in the simple dPIE model and the wide distribution 
of multiply-imaged galaxies across the cluster core reassure 
us that the absolute calibration of the total mass, Einstein 
radius and cluster substructure are reasonably accurate. 



3-4-2 Fixed Size of L* Cluster Galaxies 

We note that other galaxy-galaxy weak-lensing works have 
obtained differ ent values for r*„ t , but within the range 10 < 
r?„ t < lOOkDC (Nataraian et al. 2002:: Limousin et al 120073; 



systems. However, for two clusters, the strong lensing con- 
straints reveal the influence of additional cluster-scale dark 
matter halos: 

• A 868: The overall shape of the giant arc providing 
the constraints in this cluster is quite straight, and slightly 
curved towards the second brightest galaxy instead of the 
cluster centre. This shows that it is mostly influence by the 
presence of a secondary component, which we parametrize 
as centred on this second galaxy. 

• A 773: The large number of highly sheared arcs and 
multiple images show the presence of a dual component at 
the centre, as well as a third dark matter clump towards 
the east. Using the new spectroscopic redshifts of multiple 
images, we confirm the result of Sm05, who included 3 com- 
ponents in the lens model to explain the overall weak-lensing 
signal. 

We compare the magnification factors p computed from 
the the lens models with the photometry of each image in 
a given multiple system (last columns of Table . We find 
only a marginal agreement between the measured photom- 
etry and these magnifications, which is certainly due to the 
fact that the magnification factors are derived only at a sin- 
gle position, whereas many of the arcs are fairly extended, 
and also due to the surface brightness limits when measuring 
the photometry of the faintest images. The majority of the 
multiple images are strongly magnified, with typically /i ~ 2 
mags (or 5 — 10 on linear scale). Errors on the magnification 
reach very large values (Afj, > 1 mag) for multiple images 
in the vicinity of the critical line (/x J> 4 mags, or J> 30 on 
linear scale). 

3.4 Systematic Errors 

In this section we discuss four systematic uncertainties in 
our mass models: inability of parameterized models to fit 
the multiple-image constraints in A 1835 ( H3.4.1|l . the im- 
pact of fixing the size of L* galaxies on the modeling results 
( 1)3.4.20 . reliability of identification of substructures f £l3.4.3|) . 
and choice of parameterization for cluster-scale dark matter 
halos ( jjH3)l . 

3-4-1 Quality of Fit for A 1835 

We tried to improve the quality of the fit to the multiple- 
image constraints in A 1835 (cj ~ 3") by adding a secondary 
dPIE cluster-scale halo on the next brightest peaks seen on 
the K-h&nA light map (see H3.4.3p . as well as allowing the 
positions of the main cluster-scale halo (previously with its 
centre fixed on that of the BCG) to be free parameters. 
This did not yield any significant improvements in the qual- 
ity of the fit. The high value of <7; is due to the inabil- 
ity of the parameterized model to reproduce simultaneously 
the multiple-image systems located to the North and to the 
South of the BCG. It therefore appears that the main limi- 
tation of the current model for this cluster is the assumption 
of elliptical symmetry in the mass components. 

To explo re this, we used the ne w, more flexible method 
presented by ljullo fc Kneibl (|2009h that employs a multi- 
scale adaptive grid to refine a parametric model that is con- 
strained by a large number of strong lensing constraints - 



lHalkola et ail [2007). To quantify how much our choice of 
r *ut — 45kpc r ^13. 2 . 2[) affects mass model results, we there- 
fore ran the following test: we included rJut as a free pa- 
rameter for the 5 models having the largest number of con- 
straints (i.e. A 611, A 773, A 1413, A 1835 and Z2701), al- 
lowing it to vary in the range [10 — 100]kpc. The best fit 
parameters and quality of fit (judged by <7j) of these 5 mod- 
els are all consistent with the fiducial models in Table [6] 
specifically, the best-fit values of r* ut are all consistent with 
r cut — 45kpc within their respective statistical uncertainties. 
The degeneracy between r* ut and a* also obey the expected 
r cut cr* 2 =Const. relation - i.e. a constant amount of mass is 
assigned to galaxy-size haloes (within 10%) regardless of the 
specific values of these two parameters. This folds through to 
an uncertainty of just 15% on the cluster substructure frac- 
tions, which is negligible in comparison with the statistical 
errors. 



3-4-3 Identification of Massive Substructures 

The necessity or not to add a second or third cluster-scale 
mass component to reproduce the observed multiple-image 
systems has been determined by the number and location 
of the strong-lensing constraints. The presence of massive 
substructures has direct implications on the substructure 
measurements discussed in i}4] We therefore test the fiducial 
mass models by comparing quantitatively the distribution 
of if-band light in the cluster cores with the mass modeling 
results. 

We constructed A"-band luminosity density maps from 
the cluster galaxy catalogues discussed in t|2.2.3l with a 
Gaussian smoothing scale of FWHM = 20" (correspond- 
ing to ~ 70kpc at z = 0.2), which is the smallest scale over 
which we measured the influence of a secondary clump on 
the strong lensing images. The location of the highest peak 
(^K,max) in the JsT-band light maps always coincides (within 
1") with the centre of the cluster-scale mass distribution. 
We searched for secondary peaks in the K-h&nd light by 
adopting a threshold L > 0.5I/K,max and look for the lo- 
cal maxima within these regions (Figure white crosses). 
When comparing their location with the cluster-scale dark 
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Figure 3. if-band cluster luminosity density maps (gray scale) smoothed with a FWHM=20"gaussian. Red contours show the recon- 
structed dark matter distribution with a constant logarithmic scale. The white contours show the selection of X-band light peaks as 
local maxima (crosses, see i|3.4.3l for details). 



matter clumps, we find that every local maximum is coin- 
cident with a cluster-scale mass component in our models 
within the modelling errors (1"), except for the secondary 
peak in A 521 (located ~ 1' South of the BCG). It therefore 
appears that we are unable to detect this structure because 
the strong-lensing constraints lie exclusively North of the 
BCG. The presence of substr ucture in A 521 is co n sisten t 
with the dynamical study of iMaurogordato et al.l (2000) , 
who identified many large-scale components in A 521, and 
suggested that it is a highly disturbed merging cluster. Our 
inability to detect any of these structures may explain why 
A 521 is an outlier in the / su b — Amu relation discussed in 
34^21 

3.4.4 Parameterization of Dark Halos 

To examine the amplitude of any systematic error in clus- 
ter mass measurement arising from choice of dPIE or NFW 
mass profiles, we show the projected mass profiles Msl(-R) in 
Fig. g] from both dPIE and NFW models. The statistical er- 
ror on Msl(R) is typically 0.05 to 0.2dex, depending on the 
number of lensing constraints. We find a general agreement 
between the Msl(R) profiles obtained from the best-fitting 
dPIE and NFW models, and use the average difference be- 
tween the two as an estimate of the systematic error made 
by assuming a specific cluster-scale profile - this systematic 
error is comparable with the statistical errors. 



4 THE MASS AND STRUCTURE OF 
CLUSTER CORES 

In this section we use the strong lens models of the full 
sample of 20 clusters to construct a statistical sample of 
measurements of the mass and structure of cluster cores, for 
comparison with theoretical predictions. Despite the modest 
statistical significance achieved in these comparisons, due to 



the sample size, this is so far the largest sample of strong 
lensing clusters analyzed in a uniform manner, and repre- 
sents a four-fold increase on the strong-lensing clusters stud- 
ied by Sm05. Most importantly, these results provide a sta- 
tistical context within which to view results from detailed 
analysis of spectacular individual cluster lenses. 

We begin by describing integrated cluster masses and 
related quantities such as Einstein radii in i )4.1l and then 
concentrate on structural quantities such as substructure 
fractions and cluster ellipticities in i)4.2l 

4.1 Cluster Masses and Einstein Radii 

4-1.1 Einstein radii 

We first measure the effective Einstein radius Be of each 
cluster at zg = 2. Be is defined as the angular radius 9 
from the centre of the cluster at which the average con- 
vergence k{9 < 0e) = 1 (|Broadhurst fc Barkanall2008l ). We 
choose z s — 2 because this is the typical redshift of the spec- 
troscopically confirmed multiply- imaged galaxies (Table HJ. 
Values of 9e range between ~ 4" and 47", are listed in 
Table and the distribution is plotted in Fig. [5] The 9e 
distribution is best-fitted by a log-normal distribution with 
(log 10 9e) ~ 1.16 ± 0.28, where 9e is measured in arcsec- 
onds. Two clusters (A 1689 a nd Abell 1703) are in common 
between this study and the iBroadhurst fc Barkanal (120081) 
sample of clusters with large Einstein radii. These two clus- 
ters are located at 2a and 1.5cr above (log 10 9e), suggesting 
that ~ 2.3% and ~ 6.7% of clusters found in larger sample 
will have Be at least as large as those of A 1689 and A 1703 
respectively. 

Comparison of observed Einstein radii with theoretical 
predictions from numerical simulations is problematic be- 
cause the simulations require both sufficiently large volume 
to contain a large sample of clusters as massive as observed 
systems, and sufficient numerical resolution to allow 9e to 
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Figure 4. Integrated mass profiles assuming a dPIE (solid grey) or a NFW (hatched red) profile for the dark-matter distribution. The 
arrows mark the distance of the multiple images used as constraints in each cluster, and the vertical dashed corresponds to the radius of 
250kpc used to derive the enclosed masses. 



be measured reliably from the simulated data. Even mod- 
ern simulations such as the Millennium Simulation are un- 
able to satisfy both requirements, the main shortcoming be- 
ing the simulation volume. Nevertheless, it is interesting to 
make a comparison. We therefore take the best-fit mass and 
concentration of the ten NFW models in Table [5] and con- 
volve the measured concent rations with the predicted con - 
centration distributions in iBroadhurst fc Barkanal l|200&t l. 
We then calculate the predicted 8e distributions and sum 
them to produce the red dashed curve in Fig. [5] This dis- 
tribution is much broader than the observed distribution, 
mainly because the virial masses of the NFW models are 



poorly constrained by strong-lensing constraints alone. In 
spite of this, we find that the predicted distribution peaks 
at 8e — 5", a factor of 2 lower than the observed distribu- 
tion. This difference may be caused by an important differ- 
ence between the simulations and observations, namely the 
presence of baryons in observed universe and the absence of 
them from the simulations. However it is interesting to note 
tha t in a recent weak-l ensing study of similarly sized sam- 
ple, [OfebeelaD (|2009t) found that the normalization of the 
mass-concentration of observed clusters is roughly a factor 
of 2 higher than predicted from numerical simulations. The 



LoCuSS: Strong-lensing Clusters 17 



10 



e E (z=2) 



100 



Figure 5. The observed distribution of effective Einstein radii for 
a source at z s = 2.0 (solid black histogram), and best-fit lognor- 
mal distribution (black dotted line). We also plot the predicted 
distribution based on the NFW model parameters (Table |6j, and 
the predicted distributions of halo properties in strong-lensing 
selected simulated clusters - see ^4.1 . II for more details. 
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Figure 6. Probability distribution of lensing masses (measured 
in a projected radius of 250 kpc) for the Sm05 sample (dashed 
histogram) and the current sample (solid histogram). 



physical origin of both differences, if confirmed by larger 
samples, may be similar. 

4-1.2 Projected Mass Measurements 

We show projected mass maps of the fO new strong-lensing 
clusters in Fig. and integrate these maps and their equiv- 
alents for the extended sample to measure the projected 
mass of the cluster cores Msl(R < 250kpc) - see Table [7] 
The mass measurements for the extended sample were based 
on the published dPIE models of these clusters listed in Ta- 
bles [2] & IB H adjusted to the cosmology used in this paper 
where relevant. The aperture of R < 250kpc is chosen be- 
cause this ensures that the region within we measure mass 
lies within the HST field of view for all clusters; it is also 
~ 2x the largest observed Einstein radius within this sam- 
ple. 

We compare the Msl{R < 250kpc) distribution of the 
current sample with the 10 clusters studied by Sm05 in Fig|Bl 
and fit a log-normal distribution to both samples, obtaining: 



(log 10 (M SL )) = 14.29±0.19 and 14.27±0.17 for our sample 
and Sm05 respectively. This confirms that both studies are 
probing clusters of the same mass, which is consistent with 
close match in the range of X-ray luminosity probed: 4.3 x 
10 44 < L x < 22.8 x 10 44 ergs" 1 in the case of Sm05, and 
3.5 x 10 44 < L x < 22.8 x 10 44 ergs _1 here (Table©. 

We also plot in Fig. [7] the relationship between 6e and 
M(< 250kpc), revealing a clear positive linear correlation 
between the two quantities, as is expected from the proper- 
ties of analytic descriptions of dark matter density profiles 
from simulations. The best-fit relation for the whole sam- 
ple is 9 E = (-14.7 ± 4.7) + (13.8 ± 2.5) M SL ; the best-fit 
relations for cool-core and non-cool core sub-samples are 
both statistically indistinguishable from this relation. How- 
ever, we find that disturbed clusters tend to lie below the 
best-fit relation for the full sample - indeed, we obtain best- 
fit relations of: 9 E = (-8.6 ± 4.5) + (9.75 ± 2.3) Msl and 
E = (-14.7 ± 5.0) + (14.0 ± 2.2) Msl for disturbed and 
undisturbed clusters, respectively. We interpret this as im- 
plying that disturbed clusters have flatter density profiles 
than undisturbed clusters of comparable mass. This is likely 
due to a combination of (i) the cluster-cluster mergers that 
are likely to both soften the density profile, and also to 
cause the observed disturbance (i.e. cause the offset between 
the X-ray and optical centers of the clusters), and (ii) dis- 
turbed clusters tending to have their merger axis preferen- 
tially aligned in the plane of the sky, thus making the distur- 
bance possible to measure. This latter point is of particular 
importance with respect to A 1689, because this cluster is 
classified as undisturbed despite t here being strong evide nce 
for it being a line-of-sight merger l|Limousin et al.ll2007bl) . If 
this cluster were viewed side-on, then it would probably lie 
below the relation for undisturbed clusters, and be classified 
as disturbed. 



4-1.3 Comparison of Lensing and X-ray Mass 
Measurements 

We further investigate the integrated mass of the cluster 
cores within R < 250kpc by comparing the Msl with masses 
measured within the same aperture from the X-ray models 
discussed in t )2.4l We used least-squared minimization, tak- 
ing account of errors in both quantities to find the best fit pa- 
rameters of the relation: Msl = B Mx- For the full sample, 
we obtain a best-fit of logio(B) — 0.13 ±0.04, i.e. a mean X- 
ray/lensing mass discrepancy of a factor of Msl/Mx = 1.3 
at ~ 3ct significance. 

We then fit the same relations for the cool core/non-cool 
core as well as disturbed/undisturbed cluster subsamples, 
and report the best fit values in Table [8] We find that the 
disturbed/undisturbed sub-samples are statistically undis- 
tinguishable from both each other and from the entire sam- 
ple. However, we find that the normalization, logio(B), of 
the relation for non-cool-core clusters is 2a higher than for 
cool-core clusters. These differences are illustrated in Fig. [H] 

To test the link between the agreement between lens- 
ing and X-ray mass measurements and the mass estimates 
and the structure of the cluster cores in more detail, we 
plot Msl/Mx against a, / sub ( fl4~2~Tj) and X-ray/BCG off- 
set in Fig. [9] There is considerable scatter in all three pan- 
els, however we try to fit a straight-line relation to the 
data in each case. We obtain a positive correlation be- 
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Table 7. The mass and structure of the cluster cores inferred from the strong-lensing mass models, plus related quantities from the 
near-infrared cluster galaxy catalogs. 



Cluster 


N a 

mult 


Radial? 6 




e d 
e 2D 


e E ( z = 2) 

(arcsec) 


/sub 


^K,BCG/^K,tot 


Ami2 
(mags) 


A 521 


1 


No 


0.61±0.33 


0.66 


3.6 ±0.8 


0.13 ±0.04 


0.22 ± 0.03 


0.05 


A611 


1 


Yes 


1.76±0.33 


0.36 


21.0 ± 1.3 


0.10 ±0.01 


0.44 ± 0.06 


2.16 


A 773 


3 


Yes 


3.01±0.58 


0.385 


30.1 ± 1.2 


0.78 ±0.03 


0.13 ± 0.02 


0.13 


A 868 


1 


No 


1.97±1.11 


0.06 


14.2 ± 5.6 


0.26 ±0.12 


0.25 ± 0.04 


0.81 


Z2701 


2 


No 


1.74±0.14 


0.34 


9.0 ±0.5 


0.04 ±0.02 


0.45 ± 0.06 


2.33 


A 1413 


4 


No 


1.71±0.20 


0.64 


11.9 ±0.5 


0.07 ±0.01 


0.32 ± 0.05 


1.80 


A 1835 


7 


Yes 


2.83±0.41 


0.49 


30.5 ±0.5 


0.13 ±0.01 


0.26 ± 0.09 


1.57 


A 2204 


1 


Yes 


2.29±0.50 


0.27 


23.9 ± 2.2 


0.25 ±0.10 


0.13 ± 0.02 


0.14 


RX 11 720 


I 


Yes 


1.18±0.59 


0.59 


7.0 ±0.5 


0.10 ±0.05 


0.47 ±0.07 


1.60 


RXJ2129 


1 


No 


1.37±0.37 


0.56 


9.0 ± 1.4 


0.15 ±0.06 


0.40 ±0.06 


1.26 






Extended sample - 


see Table [2] for mass 


model references 




A68 


6 


No 


2.16±0.23 


0.23 


7.5±0.5 


0.33±0.04 


0.32±0.04 


1.40 


A383 


3 


Yes 


1.87±0.26 


0.22 


10.4±2.6 


0.02±0.01 


0.46±0.06 


1.90 


A963 


2 


No 


1.74±0.44 


0.355 


7.5±0.5 


0.13±0.07 


0.38±0.05 


1.26 


A1201 


1 


No 


0.80±0.33 


0.66 


1.5±0.15 


0.02±0.01 


0.64±0.09 


2.54 


A2218 


6 


Yes 


3.00±0.24 


0.23 


18.3±0.5 


0.54±0.01 


0.18±0.03 


0.46 


A2219 


6 


Yes 


2.33±0.23 


0.41 


15.6±0.6 


0.57± 0.04 


0.22±0.03 


0.75 


A2390 


1 


No 


1.99±0.07 


0.14 


17.5±0.5 


0.03±0.03 


0.30±0.04 


1.53 


A2667 


3 


No 


2.41±0.07 


0.39 


13.0±0.6 


0.14±0.03 


0.17±0.02 


0.87 


A1689 


34 


Yes 


4.53±0.13 


0.22 


47.1±0.5 


0.22±0.03 


0.24±0.02 


0.68 


A 1703 


16 


Yes 


2.98±0.09 


0.39 


36.8±1.5 


0.15±0.02 


0.52±0.07 


1.57 



a Number of strongly- lensed galaxies used to constrain the mass model. b Whether the strong-lensing constraints include a radial image 
pair. c Projected mass within R < 250kpc in units of 10 14 Mq. d Ellipticity of mass distribution in the cluster core inferred from the 

projected mass maps. 



Table 8. Best fit parameters for the relationship between 
M S h(R < 250kpc) and M X (R < 250kpc) discussed in 34.1.31 



Sample 


N a . 

cms 


Normalization 


Scatter 






logic (B) 


a M (dex) 


All 6 


18 


0.13 ±0.04 


0.17 


Cool core 


8 


0.05 ±0.05 


0.19 


Non cool core 


10 


0.17 ±0.04 


0.13 


Undisturbed 


12 


0.15 ±0.05 


0.20 


Disturbed 


6 


0.12 ±0.04 


0.15 



a The number of clusters in each sub-sample are taken from Ta- 
ble \E\ b Chandra data are not available for two of the clusters - 
A 868 and A 1703. 



tween Msl/Mx and all of a, / su b and offset, as shown by 
the dashed lines in Fig. [5] however the M$l/Mx vs / su b 
is the only one of the three with a slope that is statisti- 
cally distinguishable from flat, with a best-fit of Msl/Mx = 
(1.8 ± 0.3) + (0.63 ± 0.35);opi (/ sub ) - i.e. a ~ 2a detection 
of a dependence of Msl/Mx on / au b. 

We interpret this as evidence that the assumption of 
hydrostatic equilibrium required by the X-ray mass esti- 
mates is less reliable in clusters with larger substructure 
fractions. Specifically, the merging activity signaled by larger 
substructure fractions may be adding non-hydrostatic pres- 
sure support to the ICM through b ulk motions of gas in a 
manner similar to that identified bv lRasia et ail (|2006l) and 
iNagai et al.l i|2007l ) at larger radii in simulations. 



4.2 Cluster Substructure, BCG Dominance and 
Ellipticity 

4-2.1 Substructure Fractions 



Following Sm05 and ISmith fc Tavlorl ([2008), we define the 
substructure fraction / au b as the fraction of the mass within 
a radius R that is not contained in the BCG and the cluster- 
scale dark matter halo centred on the BCG. This allows 
straightforward calculation of / su b from the mass models, 
and provides an estimate of the fraction of mass that resides 
in galaxy- and group-scale halos within the parent cluster. 
We measure / su b within the same aperture as for Msl, i.e. 
R < 250kpc. At these scales / su b is stable to small changes 
in the aperture choice; this stability extends down to R > 
lOOkpc in most cases. The measurements of / sub are listed 
in Table El 



Sjnith^_Tavtor ()2008[ ) compared the predicted (from 

iTavlor fc Babulll2005l ) and observed (from Sm05) / sub distri- 
butions, finding a possible excess of clusters with the highest 
values of / su b- They speculated that, if confirmed, these ex- 
cesses may be due to X-ray selected clusters comprising an 
excess of cool core and merging clusters with respect a purely 
mass-selected sample (the syn thetic clusters were selected 
on mass from ITavlor fc Babull 's model). We return to this 
distribution in Fig.[l0]with our enlarged sample of 20 strong- 
lensing clusters (Sm05 comprises 5 strong- and 5 weak- 
lensing constrained cluster cores). As seen in Fig. 1101 the 
observed / su b distribution remains noisey with 20 clusters . 
We repeat the Kolmogorov Smirnov test o f lSmith fc Tavlorl . 
and again obtain an inconclusive result, with the null hy- 
pothesis that the observed and sythetic samples are drawn 
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Figure 9. Residuals of the Mq^/Mx relation against the central X-ray slope a (left panel), the strong-lensing derived substructure 
fraction (middle panel) or the offset between the X-ray centre and the BCG (right panel). The red dashed line shows the best fitting 
linear relation for each diagram. Symbols are identical to Fig. [8] 




Figure 10. Observed substructure fraction (solid curve), as determi ned from the 20 s trong -lensing models in the extended sample, 
compared with the expected distribution (dashed curve) presented in ISmith fc Ta ylor (200q). Left diagram: cumulative distribution. 
Right diagram: differential distribution. 



from the same underlying distribution disfavoured at 80% 
(i.e. 1.3a) confidence. 



There is extensive discussion in the literature on the is- 
sue of whether clus ter-cluster mergers are capable of des troy- 
ing cool cores (e.g. lPoole et al.ll2008l ; iBurns et al.ll200St) . We 
have previously addressed this question in directly through 
compar isons of a and X-ray/BCG offsets in lSanderson et all 
(2009a]). Here, we tackle it more directly, by comparing 
a with / su b in Fig. 1111 We fit a relation of the form 
log 10 (/ su b) = A + B a and obtain best fit parameters of 
A=0.09±0.12, and B=1.15±0.22. We therefore find a posi- 
tive correlation at 5a significance between / su b and a, albeit 
with large scatter. Clusters that host a cool core therefore 
also have less substructure in their cluster core mass distri- 
bution, and vice versa. This supports the connection found 
in H4.1.3I between the differing values of Msl/Mx for cool 
core and non-cool core clusters, and the positive correlation 
between the mass ratio and / su b- This result is also strongly 
suggestive that the cluster-cluster merger activity associated 
with larger values of / su b plays a role in destroying cool 
cores. However the appreciable scatter on this relationship 
also suggests that the physics is more complicated than a 
simple one-to-one relationship. 



4.2.2 BCG Dominance 



The dominance of the BCG over the total K-h&nA lumi- 
nosity of cluster galaxies within cluster cores was investi- 
gated by Sm05, who found a roughly monotonic relation- 
ship between Lk,bcg / Lk,tot and / su b in the sense that 
more dominant BCGs live in clusters with lower substruc- 
ture fractio ns. This has been in vestigated in more detail 
recently by ISmith et ail (l2009bl ) who showed that ~ 8% 
of 10 15 Mq clusters contain a BCG at least two mag- 
nitudes brighter than the second ranked cluster galaxy 
(Amu > 2). We build on these results here by compar- 
ing substructure in the total mass distribution in the clus- 
ter cores with the dominance of the BCGs for our enlarged 
sample of 20 strong-lensing clusters. We plot / su b versus 
both the fraction of the cluster i\-band luminosity ema- 
nating from the BCG and the magnitude gap between the 
BCG and the second ranked cluster galaxy in Fig. 1121 with 
both measured within the same R < 250kpc region as 
/ su b. Both plots show an obvious correlation, confirming 
and extended the results in Sm05. We obtain best-fit re- 
lations of: logio(f su b) = —2.80 L K ,BCG/£K,tot + 0.22 and 
logio(fsub) = — 0.53Ami2 — 0.07, both relations having a 
dispersion of a ~ 0.3dex. 
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Figure 12. Substructure fraction versus indicators of BCG dominance: fraction of the cluster core A"-band luminosity that emanates 
from the BCG (left) and the magnitude gap between the BCG and the second ranked cluster galaxy (right). Filled symbols are the new 
clusters, and open symbols are from the extended sample. 
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Figure 7. Effective Einstein radius (assuming a z = 2 source) vs 
strong- lensing mass (measured within 250 kpc). The lines show 
the general trends when least-squared fitting a linear relation for 
all clusters (black solid line), undisturbed clusters (blue dashed 
line) and disturbed clusters (red dotted lines). Red diamonds de- 
note disturbed clusters and blue triangles undisturbed clusters. 
Filled symbols are the new strong lensing models, and open sym- 
bols are taken from the extended sample. 



4-2.3 Cluster ellipticity 

We have measured the 2D ellipticity and orientation of each 
cluster core mass distribution (R < 250kpc) by fitting el- 
liptical mass contours with the IRAF routine ellipse to 
the projected mass maps based on the mass models (see 
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Figure 8. Relationship between Mx and Msli both measured 
within R < 250kpc. The lines show the general trends when least- 
squared fitting a linear relation for all clusters (black solid line), 
cool-core clusters (blue dashed line) and non-cool core clusters 
(red dotted line). Blue triangles denote cool-core clusters and 
red diamonds non-cool core clusters. Filled symbols are the new 
strong lensing models, and open symbols are taken from the ex- 
tended sample. 



also iRichard et ai1 f2009). We kept the centre of the ellipses 
fixed at the peak of the mass map and let the ellipticity 
e(a) and position angle 4>2T>{a) vary as free parameters with 
the semi major axis a. The resulting ellipticities are listed 
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Figure 11. Relationship between substructure fraction, / su b, and 
a the slope of the gas density profile at 0.04, r^oo, as a measure 
of the strength of cooling in the cluster cores core. Symbols are 
identical to Fig. [H] The dashed line shows the best-fit relation to 
the data. 



in Table [71 and the distribution is plotted in Fig. 1131 We 
also show in this fi gure the predicted distributions from 
lOeuri fc Blandfordl j2008h . both for the underlying mass- 
selected cluste r population, and fo r clusters with the largest 
Einstein radii. lOguri fc Blandfordl interpreted the difference 
between these two distributions, with the latter peaking at 
smaller ellipticities than the former, as implying that the 
major axis the 3D mass distributions of strong-lensing clus- 
ters are more likely to be closely aligned with our line of sight 
through the clusters. Strong-lensing clusters are therefore 
expected to be rounder on the sky than non-strong-lensing 
clusters. 

Our observed distribution agrees better with the pre- 
dicted distribution of large Einstein radii clusters than with 
the mass-selected population, however we detect a peak 
at e2D — 0.65. This peak contains A 521 and A 1201, the 
two clusters in our sample with the smallest Einstein radii, 
demonstrating that clusters with small Einstein radii can 
still be strong lensing clusters if background galaxies are 
fortuitously aligned. We therefore exclude these two clusters 
from the observed distribution and fit a normal distribution 
to the other 18 clusters, obtaining: (e2o) = 0.34 ± 0.14 



5 SUMMARY AND DISCUSSION 

We have presented a strong lensing analysis of 20 massive 
galaxy clusters from the X-ray selected LoCuSS sample, all 
of which contain at least one spectroscopically-confirmed 
multiply-imaged background galaxy. Ten of the sample are 
new strong-lensing clusters for which we present the first 
detailed parametric lens models based on new spectroscopic 
redshift measurements at the Keck I telescope. All of the 




e2D 

Figure 13. Observed distribution of 2d ellipticities (red solid 
histogram) compa r ed wi th predictions from simulations by 
lOguri fc Blandfordl l|2008h . for an unbiased cluster population 
(black dashed histogram) and a cluster population producing 
the largest Einstein radii (black dotted histogram). The observed 
peak at e2o ~ 0.65 is dominated by A 1201 and A 521, the 2 clus- 
ters with the smallest Einstein radii (red dot-dashed histogram). 
The red solid curve shows the gaussian fit to the probability dis- 
tribution. 



clusters are well described by a parametric mass distribu- 
tion containing one or several cluster-scale dark-matter halos 
centred on the brightest peaks of the cluster K-band light. 
We used our parametric models to compute maps of the 
mass distribution in the cluster cores, and thus to measure 
the projected mass and fraction of that mass associated with 
substructures within R ^ 250kpc. We have compared these 
measurements with observations of the baryons within the 
clusters, employing X-ray and K-band data to probe the in- 
tracluster medium and evolved stars respectively. Our main 
results are summarized below: 

(i) The Einstein radius for a typical source redshift of z — 2 
spans 4" ^ 0e(z — 2) ^ 47" and is lognormally distributed, 
peaking at (log 10 9e) = 1.16 ± 0.28, where the uncertainty 
is width of the distribution. Famous clusters with a large 
Einstein Radius, A 1689 and A 1703, are thus "outliers" at 
2a and 1.5cr above the mean respectively. We also compare 
the observed distribution with that predicted from the Mil- 
lennium Simulation, and find that the peak of the predicted 
distribution lies at ~ 5", i.e. ~ 1.7c below that of the ob- 
served distribution. 

(ii) 6 E (z = 2) is correlated with M S l(R < 250kpc), the 
projected mass of the cluster core obtained from the strong 
lens models. We find that "disturbed" clusters, i.e. those 
with an offset between the centroid of their X-ray emission 
and the optical centroid of their BCG of > O.Olrsoo, typically 
lie below the best-fit 9e — Msl relation. We interpret this 
as arising from a combination of (i) cluster-cluster mergers 
(assumed to be responsible for the X-ray/BCG offset) acting 
to soften the cluster density profile, and thus reduce 0e, and 
(ii) an orientation effect, i.e. that the major axis of disturbed 
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clusters tends to be closer to orthogonal to the line of sight 
through the cluster, whilst undisturbed clusters tend to have 
their major axis parallel with the line of sight. 

(iii) The ratio of strong-lensing- and X-ray-based projected 
cluster mass measurements within R < 250kpc is measured 
to be Msl/Mx = 1.3, discrepant with unity at 3cr. This 
X-ray/lensing mass discrepancy depends on the structure of 
the cluster core - we detect a positive correlation between 
Msl/Mx and the fraction of cluster mass associated with 
substructures within R < 250kpc, / su b, at 2a significance. 
We interpret this as evidence that the cluster-cluster merger 
activity associated with cluster substructure is responsible 
for departures from hydrostatic equilibrium. 

(iv) The substructure fraction, / su b, is also correlated with 
a, the slope of the logarithmic gas density profile at 0.04rsoo, 
in the sense that clusters with steeper (more nagative a) gas 
density profiles have smaller values of / su b, and vice versa. 
The gas density profile slope is used as an indicator of the 
strength of cooling in the intracluster medium, the steepest 
slopes (a < 0.85) being identified as "cool core clusters". 
This direct empirical relationship between a and / su b, im- 
plies a connection between the cluster-cluster mergers and 
the strength of cooling in cluster cores. 

(v) We also find a strong correlation between / su b and the 
dominance of the BCG in the sense that clusters with more 
substructure have less dominant BCGs. This suggests that 
measures of BCG dominance, including the luminosity gap 
statistic (difference between the magnitudes of the first and 
second ranked galaxies), may be a useful probe of cluster 
substructure. 

(vi) The distribution of cluster core ellipticities, measured 
from the mass maps that are in turn computed from the 
parametric lens models, are consistent with the predicted 
ellipticity distribution of strong-lensing clusters. 

This is the largest published sample of strong-lensing 
clusters to date - a 4x increase in sample size from 
ISmith et al.1 (|2005l) . Overall our empirical results are con- 
sistent with those of ISmith et al.l in that a straightforward 
interpretation is that cluster-cluster mergers play a promi- 
nent role in shaping the observed properties of cluster cores. 
This is most striking in the X-ray/lensing comparisons, i.e. 
significant detection of a 30% X-ray/lensing mass discrep- 
ancy, and a dependence of this discrepancy on structure of 
the cluster core. This is underlined by the correlation be- 
tween / su b and a - implying that cluster-cluster mergers 
both cause departures from hydrostatic equilibrium and play 
a role in moderating the cooling of gas in cluster cores. This 
latter point rema ins controversial (e.g. iPoole et al.1 120081 ; 
iBurns et al.ll2008h . and we caution that the distribution of 
formation epochs, and other physical processes such as pre- 
heating may ultimately modify our simplistic interpretation. 

We also note that the clusters with the largest Einstein 
radii in our sample (A 1689 and A 1703) are <J 2a above the 
mean of the best-fit log-normal distribution, implying that 
~ 2 — 7% of larger samples will contain comparable clusters. 
We also speculate that the ~ 1.7<r discrepancy between the 
peak of the observed and predicted Einstein radius distribu- 
tions may be in part attributable to the presence of baryons 
in the observed universe, in contrast to the dark matter only 
simulations on which the prediction is based. Certainly, the 
discrepancy is likely over-estimated because of the current 



computational limitations on numerical simulations, namely 
that even the most advanced simulations such as the Millen- 
nium simulation do not embrace sufficient volume at suffi- 
cient numerical resolution to contain a representative sample 
of the most massive strong-lensing clusters. In summary, we 
find no compelling evidence from our statistical analysis of 
a sample of strong lensing clusters to support recent claims 
that clusters with lar ge Einstein radii present a ch allenge to 
the CDM paradigm l|Broadhurst fc Barkanall2008l ). 

This article has concentrated on clusters observable 
from Mauna Kea with the Keck-I telescope. In the future we 
will expand this sample to include Southern clusters that we 
have also observed with HST and followed up spectroscop- 
ically with VLT and Gemini-S (P. May et al., in prepara- 
tion) . Similar st udies of clusters at z > 0.3 from the M Assive 
Cluster Survey labeling et al.ll200ll . 120071 1 will allow evolu- 
tion in the proper ties of cluster cores to be probed (e.g. 
ISmith et al.ll2009al Richard et al. in preparation). Within 
LoCuSS, future work will include combining these strong 
lens models with our recentl y published weak-le nsing anal- 
ysis of Subaru observations l|Okabe et al.l [2009h - this will 
allow a detailed investigation of the structure of the clus- 
ter mass distributions across a wide range of physical scales. 
The well-constrained lens models presented here are also 
well-suited to be utilized in grav itational telescope searches 
for very high red shift galaxie s (IMaizv et al.l 120091) . f ollow- 
ing, f o r example. iKneib et al.l (|2004h : iRichard et all (|2006l. 
2008) ; iBouwens et al.l (|2009h . Finally, we have presented op- 
tical spectra of lensed background galaxies at z ~ 1 — 4, 
which are magnified by 1 — 4 magnitudes from strong lens- 
ing. These sources are well suited for further high resolution 
spectr oscopic follow-up , such as near-infrared IFU observa- 
tions l|Swinbank et al, 



p, suc n 
2009). 
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APPENDIX A: SPECTRA OF MULTIPLE 
IMAGES 
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Figure Al. Extracted spectra of multiple images from the current sample. We mark the proeminent spectral features used to derive 
the rcdshift. 
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Table Bl. Best-fit parameters of the mass models. For each mass component, we give the centre, ellipticity, orientation, core and cut 
radii, as well as central velocity dispersion of the dPIE profile. The following column gives the image plane RMS of this model 
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